Collapse and revivals of the photon field in a Landau-Zener process 
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We consider the evolution of a two-level system coupled to a photon field initially in a coherent 
state, as the energy of the two-level system is linearly varied through resonance with the photon 
field. At a fixed time after the resonance, the amplitude of the photon field is found to show a 
collapse and subsequent revivals as a function of rate of energy variation. Including decay of the 
photon field, we find that the observation of such collapse and revivals is near the technological limit 
of current cavity QED experiments but should be achievable. 
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The famous problem of Landau [l[ and Zener Q con- 
cerns a two-level system whose parameters are varied so 
that an anticrossing of energy levels occurs, and provides 
the probability that the system will remain in an adia- 
batic state. There has been much work since on finding 
generalisations of this problem to many levels which can 
still be exactly solved, e.g [1, 0, 0, @, 01 • Another generali- 
sation of the Landau Zener process is to consider multiple 
occupation of single particle levels, and then to find how 
varying the single particle level energies affects the many 
particle statel8| . These constitute a large class of possible 
problems, of which wc study a particular case relevant 
to cavity QED experiments % [3, [HI, [H, [H. Cavity 



QED studies strong light-matter coupling between pho- 
tons confined in a cavity and matter (e.g. single atoms, 
or artificial atoms such as quantum dots [l^ or Joseph- 
son junctions We consider a cavity QED system in 
which the matter is driven through a Landau-Zener tran- 
sition, and find that this has dramatic consequences for 
the state of the photon field: observing at a fixed time 
after the anticrossing and changing the rate of energy 
variation, there is a collapse and revival of the photon 
field amplitude. These oscillations extend far into the 
regime of slow energy variation — the adiabatic limit 
— in which the single particle Landau-Zener probability 
shows no further change. 

The scheme we propose consists of a two-level system 
coupled to a photon mode in a cavity; the energy of 
the two-level system is varied linearly in time, passing 
through resonance with the photon energy. A related 
problem, but with a classical photon field, was consid- 
ered in Ref . [3| ; we also note that a generalisation of our 
model to many two-level systems is analogous to mod- 
els of interest describing production of molecules in cold 
atomic gases, when varying the molecular energy by a 



Feshbach resonance, e.g.[8j, [16|, [17|, [18[- If the initial cav- 
ity state were a number state then the result would be 
the Landau-Zener energy level crossing: The state with 
an excited two-level system and n — 1 photons crosses 
the state with an unexcited two-level system and n pho- 
tons; the mutual repulsion of these states is given by 
g\/n, where g is the radiation-matter coupling strength. 



Varying the two-level system energy linearly in time as 
E = Xt induces transitions between these states. The 
standard Landau-Zener result [l|, 0] is that in terms of 
the parameter z = the probability of remaining in 

the lowest energy state (rather than exciting the system) 
is P = 1 — e"^'^^". If the energy ramp is slow, or zn^ 1, 
this quickly saturates at 1. 

However, a natural state of the cavity is the coher- 
ent state [19|, leading to a far more intricate behaviour. 
The coherent state can be resolved onto a basis of num- 
ber states, and each number state evolves independently. 
The effective radiation-matter coupling strength for n 
photons is g^/n, and so the evolution of each photon 
number state is also n dependent. In particular, each 
photon number state acquires a different phase during 
the process; this phase difference continues to play an 
important role even when the sweep velocity is so slow 
that z ^ 1 and the probability for a single transition 
P = 1 — e^^^^" ~ 1. To measure the effect of these phase 
differences, the easiest quantity to measure is the ampli- 
tude of the resultant photon field, {tjj) = How 
this is to be measured depends on the particular cavity 
QED system used. Some possible methods include Ram- 
sey interferometry of a second atom to probe the state of 
the cavity [l3l |: homodyne measurement of the photons 
leaking out of the cavity by interfering it with a reference 
beam; and homodyne measurement inside the cavit y by 
detecting the interfering fields with a second atom [12j . 
The dependence of this field amplitude on the rate of 
sweeping is shown in Fig. [l] for two different values of 
the initial field amplitude. 

The photon field amplitude shows a collapse and subse- 
quent revivals arising from the phase difference between 
different number state components of the coherent state. 
The collapse occurs as, with slower driving, the phase 
difference between different number states grows, and so 
their contributions to the field amplitude interfere de- 
structively. At yet slower driving, the phase difference 
grows sufficiently that subsequent number states are back 
in phase. However, since the dependence of the phase on 
photon number is nonlinear, these revivals are imperfect, 
and eventually cease. 
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FIG. 1: Dependence of the photon field intensity on the 
Landau-Zener parameter, z — /\. Main figure, initial co- 
herent state amplitude i/; = 5; inset ^ — 10. The red dashed 
line is the approximation discussed in the text. 



Let US now discuss this behaviour more quantitatively; 
the Hamiltonian of the problem is: 



H 



At the beginning of the time evolution, the two-level 
system energy is well below the photon energy and 
the system is initialised in a coherent photon state, 
l^f,) = e-l'^l'/2^„(V'"/y^) |n,t) . Evolution under 
Eq. H]) gives \n, t) ^ A^+i \n, j) +B„+i \n + 1, j), where 
for evolution from — T to T for large times T, the 
coefficients 0, [i^l are: 



An = e" 



-mrz /2 



r(l — inz) 



'i[Tr/4+\T^ /2+2znln(,^/\T)] 



(2) 
(3) 



In the large z limit where interesting behaviour is seen 
in Fig. [1] An — and ~ 1 so we may write i?„ = 
exp(i0„) with 4'niT) = zn [ln(2;n) — 1] — 2znhi(\^T) + 
(j)Q. In this adiabatic limit, this phase is just the in- 
tegral of En{t), the instantaneous energy of the eigen- 
state with n excitations. The final many body state 
is then: |*) = e^l^l'/^ E„(V'"e*^"+i(^VV^) I" + 1, i) ■ 
The measurable field amplitude can then be written as 
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(4) 



The existence and explanation of the collapse and re- 
vivals here is related to the collapse and revival of Rabi 
oscillations 21, 22, 2^; this occurs in a model like Eq. ([1]) 
but without time- varying energies, and collapse and re- 
vivals occur as a function of time. Let us discuss the 
case IV"!^ ^ 1. As in |2l|, one can expand the phase 



difference A(/>„ = (j)n+i — 4>n near n = (where 
the amplitude of the terms in the sum peaks), giving 
A4>\.4,\2+„, = ^^[W + zm/\if - 2mV(2|V'|^). The re- 
vivals in Fig. [1] occur when A0„ is the same (modulo 
27r) for each term in Eq. (|3]); for small m this condition 
is z = 2'kN\iI)\'^, where N is an integer labelling the re- 
vival. Near such revivals, after subtracting 2'iTNm from 
A(/)|^|2_|_„, it can be written so it varies slowly with m: 



A(/)|.0|2+r„ - 27riVm = A0|^|2 + 



(z - 27riV|V'nm zm^ 



This then allows the sum in Eq. (j4|) to be replaced with 
an integral over m = n — |f/'P- Using the Gaussian ap- 
proximation to the Poisson distribution, yields: 



dm 



exp 



+ «A0|^|: 



Evaluating this Gaussian integral, and summing over val- 
ues of N for each revival gives the result: 



IV'I 



z^exp 

N=0 



-{z -2ti N\i:\'^f 
2|^|2(l + z2/|V'l^ 



(5) 

At very large z, the revivals disappear and the behaviour 
becomes complex because terms of higher order in m in 
the expansion of A0|.0|2^„ play a role when z ^ 3|V'p. 
This means revivals are not seen for N > TVmax ~ I "01/2- 
Equation ([5]) is shown in Fig. [T] by the red dashed line. 

The revivals seen in the field amplitude do not in- 
dicate a complete revival of the initial coherent state. 
This can be seen by considering the Wigner function [l^ , 
W{x,p) = / **(a; -h y)^{x - y)e^'Py dy / tt , where ^{x) = 
{x\^) is the position representation of the photon wave- 
function, and can be written in terms of Hermite poly- 
nomials Hn{x) (which correspond to number states) as: 



l'(x) = 



oo 

■E- 



Hnix). (6) 



As seen in Fig. [21 as soon as z ^ 1 the Wigner func- 
tion has many nodes, and describes a highly non-classical 
state, unlike the initial coherent state whose Wigner func- 
tion is a Gaussian [l^. We note that unlike the collapse 
and revival of Rabi oscillations 2^, |23| this non-classical 
state is not associated with entanglement between the 
two-level system and photon field [2J] at the end of the 
dynamics, although transient entanglement does occur 
during the sweep. Since the collapse and revival occur 
in the adiabatic limit the two-level system always ends 
in a pure state and so (in the absence of decay) the pho- 
ton state is a pure but non-classical state. The complete 
Wigner function can be measured by quantum state to- 
mography 2^, 2^ ; such measurements have recently been 
performed in cavity QED experiments [S^, . 
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FIG. 2: Wigner function of the final photon state for ^ = 5 at 
values of z corresponding to the first collapse and first revival. 
Inset: Glauber Q function on same scale for same parameters. 



In a real cavity QED experiment there is decoherence 
due to escape of photons out of the cavity, and decay 
of the two-level system without emitting a photon; these 
may present an obstacle to observing the revivals, and so 
we next investigate the maximum photon loss rate k and 
non-radiative decay rate 7 that can be tolerated. With a 
non-zero decay rate k, one must consider an finite dura- 
tion of the level crossing, and to maximise the signal one 
should make the duration as short as possible consistent 
with accumulating the necessary phase. A naive estimate 
of the effect of decay comes from {nt* ,jt*) where t* is the 
time for the level crossing to occur, given by Xt* ~ <?|'0li 
and (7 1 ■01 is the characteristic coupling strength for the 
coherent photon state. Decay is weak if {Kt*,jt*) ^ 1, 
which gives {K,j)/g ^ (A/g^)/|-0|. For this condition to 



be satisfied at the first revival, when g /X 



2tt\tP\' 



one requires that {K,j)/g ^ l/(27r|V'P)- Away from the 
level crossing the effects of decay are less serious; decay 
after the crossing attenuates the amplitude of the sig- 
nal, decay beforehand reduces the amplitude of the initial 
photon field amplitude. This naive estimate is sufficient 
for the effects of non-radiative decay, but not for loss of 
photons. 

To investigate the effects of decay quantitatively, we 
solve numerically the density matrix equation of motion: 



dtp 



H,p 



(7) 



where £k[/5] = {k/ 2)[^(;^ ijj p+ f/'ip''' iIj — 2'tpptp'''] describes loss 
of photons and Cj[p] = {j/2)[a^aop + pa+a^ — 2a_pa^] 
describes non-radiative decay of spins [l9|. The numer- 
ical solution of this equation, shown in Fig. [31 reveals 
a greater effect of photon loss than the naive estimate 
above. For -0 = 5, the maximum permissible decay rates 
is K/g ~ 2 X 10~^, 7/(7 ~ 2 X 10~^, compared to the naive 
estimate I/(27r|-!/'P) ^ lO'^. 

To understand this enhancement of decay, it is conve- 
nient to calculate the effect of photon decay in the adia- 
batic limit, such that the state \n — I, t) evolves to: 

(8) 
(9) 
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FIG. 3: Effect of decay on z dependence of field amplitude. 
Left panel; effect of photon leakage, k. Right panel: ef- 
fect of non-radiative decay 7 (with field amplitude plotted 
on logarithmic scale). Plotted for ?/j = 5 and a time range of 
-10" <gt< 10*. 



The adiabatic approximation is appropriate even with 
photon decay because the probability that loss of a pho- 
ton swaps between the |-|-) and |— ) subspaccs is small. 
This probability, of switching from -f) to \n — 1,—) 
is given by: P_+ = |(n - 1, -l^jn, +10^^!"-: +) 
and can be written in terms of the t of Eq. as: 



3(6'„) sin(6'„_i)^ - sin(6'„) cos(6l„_i)v'n - 1| 



16n2 



sin''(6'„ 
1 2 



< 



1 27 
16n2 16 



(10) 



where the last expression makes use of approximations for 
n ^ I. It is thus clear that for large n, the probability 
of leaving the adiabatic subspacc is small. 

In this adiabatic subspace, the density matrix equa- 
tion gives a closed set of equations for off-diagonal terms. 
An = (n — 1, -f |/5|n, and the measured field amplitude 
is (*|'!/'|*) = Y^nVnAn- The effect of the Hamiltonian 
is phase evolution of each A„, with the final phase gain 
of A„ being A</)„_i(T) as in Eq. (j4|). Including also the 
matrix elements due to photon decay, one can write: 



dAn 
dt 



- i j A„ - ^fj 



n+l 



(11) 

For the range of parameters shown in Fig. [31^ a), the re- 
sults of numerical evaluation of this equation and of the 
full problem, Eq. ([7]) cannot be distinguished by eye. 

If K = 0, the solution for A„ recovers Eq. ([4]). Al- 
ternatively, if A0„ is time independent, then for an 
initial coherent state, the time evolution can also be 
exactly solved [l^ and the field amplitude decays as 
ip(t) = -0(0)6^"*/^. It may appear surprising that the 
characteristic decay rate is k, while the apparent rate in 
Eq. ()lip is nn; the explanation is that the order n con- 
tributions from the A„ and A„+i terms cancel. However, 
with a time-dependent A0ji, the A„ and A„_|_i terms pick 
up a relative phase difference, and the cancellation of the 
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order n term fails. Thus the characteristic decay rate 
is Kn ~ during the level-crossing time, hence the 

actual requirement for small decay is ^/g <C l/(27r|?/;|^); 
for = 5 this gives k/^ <C 5 x 10~^ (cf Fig. [3]). 

Of the various realisations of cavity QED, those with 
parameters closest to those required by the constraints 
on K/g,^/g are either Josephson junctions coupled to 
stripline resonators or Rydbcrg atoms in microwave 
cavities [2^. For Rydberg atoms the values oi n/g, 7/5 
currently achievable are very close to those required, but 
the limiting factor is the short time taken for an atom to 
pass through the cavity. 

Finally, let us present a more quantitative understand- 
ing of the "enhanced decay" by comparing the perturba- 
tive solution of Eq. pT|) to the naive expectation: 



can be found in terms of a Besscl function yielding: 



«'|^(k,Vo)|* 



(12) 

This naive decay describes decay of the initial coherent 
state before the anticrossing (thus shifting the revivals 
to smaller values of z), and decay of the field amplitude 
after the anticrossing. The perturbative solution (to lead- 
ing order in k) of Eq. (|lip can be easily found by writing 
A„ = A„ exp[— iA(/)„_i(t)], and then ignoring the n de- 
pendence of A„ on the right hand side of Eq. pT|) . It 
will be convenient to write the initial conditions for A„ 

as A„ = P„V/V", where P„ = e-l'^l'|V'p("-^V(" - 1)! 
is the probability of having n — 1 photons. Writing 
(5('I'|-0|'I') = ('J'|-0|'J')naivc - (^'IV'I*) and expanding both 
Eq. and the solution to Eq. pT|) to leading order in 
K, one has: 



jigi[A0„_i(-T)-A0„_i(T)] _^ 2igi[A0„(-T)-A0„(T)] 



jj[A0„(-T)-A0„_i(T)] / gi[A0„_i-A0„]^^ 



(13) 



Assuming n 3> 1, one can simplify the expressions for the 
phase differences appearing here to give: 



dt 



(-) 

\2nJ 



2n + At2/4 ] 



(14) 



Finally, note that the integral in Eq. (11411 has a finite 
limit as r — *■ 00; this means that the naive decay that 
was subtracted fully describes the decay at long times, 
and Eq. (|14p describes only the extra contribution that 
occurs during the crossing. Taking T 00, the integral 



TTX 




4^1 ff) (15) 
?iA \2nJ ^ ' 



Using the asymptotic form for the Bcssel function, this 
expression shows that the characteristic scale of the 
extra decay is given by: (5(^|'i/'|^) oc k\iP\'^iIjz/^/\ cx 
(K/g)|'0pz'^/2. This extra term explains both the scale 
of the extra decay observed in Fig. [3ja), and also the z 
dependence of the decay visible in that figure. 

In summary we have proposed an experiment in which 
collapse and revivals of a coherent field amplitude can 
be seen as a function of varying sweep rate in a Landau- 
Zener level crossing problem. Considering a leaky cavity 
and non-radiative decay, this effect continues to survive 
as long as decay is weak enough, however the sense of 
"weak enough" , K/g < 10~^, differs from the naive expec- 
tation. As such, this proposed experiment is at the limit 
of the capability of current experiments, and so provides 
a dramatic consequence to be observed by any system far 
in the strong-coupling limit 
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